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Abstract 

Using a conventional BCS - Fermi liquid model we calculate the c-axis optical 
reflectivity of the layered high temperature cuprate superconductors by ob- 
taining the finite temperature dynamical dielectric function in a microscopic 
self-consistent gauge invariant formalism. We get good semi-quantitative 
agreement with all the existing experimental data by using the measured 
normal state dc resistivities as the input parameters in obtaining the c-axis 
hopping amplitude and the normal state level broadening in our microscopic 
calculation. 
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The c-axis charge dynamics of high-temperature layered cuprate superconductors is un- 
usual and interesting, and has therefore attracted a great deal of well-deserved attention. In 
this Letter we present a microscopic theoretical calculation for the c-axis optical reflectivity 
of the cuprates, concentrating on the optimum doping (normally corresponding to the high- 
est critical temperature T c ) single layer (eg. LSCO, Tl-2201) cuprate systems which have 
only one CuO layer per unit cell. An impressive series Jl]-[3| of c-axis optical reflectivity 
measurements by Uchida and co-workers on single-crystal Lai-xSrxCuO^ (LSCO) systems 
shows the following salient features in the frequency dependence of the electronic reflectiv- 
ity: (1) For temperatures T > T c , i.e., in the normal phases, the c-axis optical reflectivity is 
essentially structureless down to 20 cm" 1 (the lowest frequency studied in ref. U), and the 
data resemble that of an ionic insulator; (2) in the superconducting phase, for T < T c , the 
reflectivity develops a well-defined sharp plasma edge with the effective plasma frequency, 
ujpi (defined as the frequency where the real part of the dielectric function goes through a 
zero, Ree(cOpz) = 0), being below the superconducting energy gap A (in fact, u p i ~ 50cm -1 
in optimally doped LSCO) - thus the superconducting state exhibits a classic "metallic" op- 
tical reflectivity behavior in contrast to the normal state "insulating" behavior; (3) the very 
low c-axis plasma frequency (below the energy gap) associated with the observed reflectivity 
plasma edge corresponds to a carrier mass substantially larger than the c-axis optical mass 
obtained in standard band structure calculations, i.e., the experimental u p i is much smaller, 
by more than an order of magnitude, than the corresponding local-density-approximation 
(LDA) band structure calculated value; (4) the superconducting state plasma edge shows 
a characteristic temperature dependence, becoming weaker (and uj p i(T) shifting to lower 
frequency with increasing temperature) as T approaches T c from below; (5) the plasma edge 
shows a strong doping dependence. It has been emphasized that these anomalies, in partic- 
ular, items (l)-(3) above, are very difficult to understand on the basis of the standard Fermi 
liquid BCS-type generic theories of high temperature superconductors. 

In the theoretical work presented here we obtain the c-axis reflectivity of the layered 
cuprate materials by calculating the appropriate long wavelength dynamical dielectric func- 
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tion e{uS) of the system. In calculating the long wavelength frequency dependent dielectric 
function we use the Nambu-Gorkov formalism Q and carry out a fully gauge invariant self- 
consistent finite temperature linear response calculation H including "impurity scattering" 
induced level broadening effects in the theory ||. Our theory, therefore, corresponds to the 
one H developed by Fertig and Das Sarma for anisotropic layered superconducting systems 
except for three important generalizations beyond ref. ||: (1) we use a finite temperature 
formalism in contrast to the T = case studied in ref. [||], (2) we use the d-wave order 
parameter symmetry for the superconducting ground state in our response calculation @ 
in contrast to the s-wave model used in ref. || (this turns out to be unimportant in un- 
derstanding reflectivity measurements), and (3) we include level broadening effects in our 
microscopic calculation in contrast to ref. || which neglected disorder. We point out that 
the d-wave symmetry makes negligible difference with respect to the s-wave symmetry in 
our long wavelength dielectric response calculation (there are substantial differences |7j at 
large wave vectors, which are not of any significance to optical reflectivity measurements). 
Inclusion of level broadening in the dynamical dielectric function is, however, crucial in un- 
derstanding the c-axis reflectivity measurements - in particular, strong disorder scattering in 
these systems, as reflected in their anomalously high c-axis dc resistivities for T > T c , gives 
rise to an overdamped plasmon in the normal state, completely suppressing any normal state 
metallic plasma edge for T > T c . A plasma edge in the superconducting state corresponding 
to the c-axis Anderson-Bogoliubov plasma mode H0] is, however, preserved because single 
particle level broadening does not affect the condensed carriers in the superconducting state. 

We assume a simple single tight binding band (of bandwidth 2t c ) along the c-axis 
and take the single-particle energy dispersion to be E{k, k z ) = h 2 k 2 /2m — t c cos(k z d), where 
k = \k a b\ is the wave vector in the ab plane (with m as the planar effective mass for charge 
dynamics in the ab plane), d is the c-axis layer separation, and k z the wave vector along 
the c-axis. Following refs. and [JF] we carry out a self-consistent (and gauge invariant) 
linear response calculation for the dynamical dielectric function e(u) treating the long range 
Coulomb interaction (in the highly anisotropic layered system) in the random phase approx- 
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imation ("bubble" diagrams) and the BCS-like short-range pairing interaction in the ladder 
approximation (which is a conserving approximation because the BCS self-energy producing 
the superconductivity is a Hartree-Fock single-loop self-energy). We include level broaden- 
ing effects both in the single-particle Green's function and in the vertex correction through 
the ladder impurity diagrams |J. Assuming the impurity scattering to be short-ranged and 
isotropic, we can parameterize the disorder strength by a single level broadening parameter 
(i.e., an imaginary part of the self-energy) 7 where r = (2 / y)~ 1 , with h = 1 throughout, is 
the transport relaxation time for the c-axis normal transport properties. Thus in principle, 
the normal state c-axis resistivity p c of the system defined 7 within our model. The calcu- 
lated microscopic e{uS) therefore depends on 7, T, A and T c (which are determined [[|,[7| by 
the strength of the pairing interaction in our BCS model) as well as on the effective single- 
particle parameters m and t c . The results depend also on a number of obvious parameters 
such as the two-dimensional carrier density n2D (which determines the density of states), 
the layer separation d in the c-direction, and the effective background dielectric constant k. 
Our philosophy at this stage is to treat 7 and t c , the two most contentious parameters in 
our theory, as unknown parameters, and take the other parameters from the best available 
experimental (and where appropriate, from band structure) results. We note that the details 
of the pairing interaction or the impurity potential do not affect our reflectivity results in 
any qualitative way. 

In Fig. |l| we show our calculated c-axis optical reflectivity to be compared directly with 
the experimental results of Tamasaku et. al. jl). We have use parameters corresponding 
to LSCO except for t c and 7 which are varied to produce different curves in Fig. |l| Our 
calculated results look very similar to the experimental observations (see, for example, Fig. 
1 in ref. []!)), and in fact our theory agrees with the measurements quantitatively if we treat t c 
and 7 as phenomenological parameters and choose t c ~ O.lA — 0.3 A and 7 > A. Typically 
we find the plasma edge u p i « 2t c — At c in the superconducting state, and provided 7 > A 
no plasma edge shows up in the normal state with the normal state reflectivity remaining 
essentially constant around 0.5 (exactly as in the experimental results of ref. JT|) down 
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to rather low frequencies (u < 0.1 A). We predict that if the normal state reflectivity 
measurement is pushed down to rather low frequencies (~ 10cm _1 or below), one would see 
a rise in the reflectivity as it would increase to unity to be consistent with there being a 
finite (albeit very low) dc conductivity along the c-axis. 



measurements |]-^| within our model one must have a very low (high) value of the c- 
axis hopping amplitude t c (level broadening parameter 7). In particular, our theory gives 
uo p i ~ 2t c in the superconducting state, which constraints t c to be in the range of 20 - 30 cm -1 
in LSCO for the theory to be quantitatively consistent with the experimental measurement 



experimental anomaly [1] of the nonexistence of any reflectivity plasma edge in the normal 
state is "explained" in our theory by having a large damping or broadening term 7 (> A) 
which enters our microscopic impurity scattering diagrams in the dynamical polarizability 
function. We have so far treated t c and 7 as phenomenological parameters, and the crucial 
issue therefore is to justify (or, at least provide a rationale for) having such small (high) 
values of t c (7) in layered high temperature superconducting systems. This, in fact, was the 
basis of the trenchant analysis of the experimental results carried out by Anderson || , 
which led him to conclude || that a standard Fermi liquid-BCS type analysis (exactly of the 
type we carry out in this work) is incapable of explaining the experimental c-axis reflectivity 
data. 

In this work we treat t c and 7 as phenomenological parameters to be fixed by normal 
state dc transport properties of the systems. We use the measured values of normal state p c 
and pab at T > T c to fix t c and 7 within the standard Fermi liquid picture. The hopping 
amplitude t c is then given simply by the formula 



The level broadening 7 is given simply by the c-axis dc resistivity through the usual relax- 
ation time approximation in the tight binding limit: 



It is clear from the results shown in Fig. [I] that to get agreement with the experimental 



in ref. QTJ] which finds a sharp reflectivity plasma edge in the 30 - 60 cm 



. The other striking 




(1) 
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e 2 n 2D dt c p c e 2 / — 5 — 

7 = 2 = ^V™*dP°P>*/ 2 - ( 2 ) 

The simple Drude-Boltzmann-Fermi liquid formulae given by Eqs. ([[]) and @) immediately 
explain why the hopping parameter t c (the level broadening parameter 7) is extremely small 
(large) - both of these anomalies arise from the anomalously large c-axis resistivity, which 
is a universal feature for the layered cuprates. We should emphasize two significant aspects 
of the phenomenological analysis giving as the all-important parameters t c and 7 (from the 
normal state dc resistivities of the system): (1) standard LDA band structure calculations 
|| give values of t c which are one to two orders of magnitude larger than those given by 
Eq. ([!]) for layered cuprates (while the two dimensional electron density and the effective 
mass in the ah plane seem to be reasonably, i.e., within a factor of two, well-described 
by LDA band calculations ||); (2) our Drude- Fermi liquid phenomenological description 
characterized by Eqs. (|1|) and (0) is completely independent of the actual normal state 
transport mechanism along the c-axis and depends only on the crucial observation that the 
normal state dc transport along the c-axis is "metallic" around optimal doping (albeit with 
an anomalously large imaginary part of the self-energy 7). 

In Table 1 we summarize our results for a number of layered cuprates, giving the values 
of the parameters t c and 7 ( obtained from experimental p c and p a b at T > T c ) used in our 
calculations to obtain u p i at T = 0, the long wavelength c-axis Anderson-Bogoliubov plasmon 
frequency in the superconducting phase. In carrying out our microscopic calculation for the 
"multilayer" cuprates (e.g. YBCO, BISCO), which contain more than one CuO layers per 
unit cell, we are concentrating here only on the low frequency c-axis plasma properties which 
arise from the intercell hopping, neglecting the mode mixing effects arising from the intracell 
hopping between the CuO layers. This is equivalent to making an effective single layer model 
for these multilayer cuprates, which should suffice for the c-axis reflectivity properties. From 
the available experimental results given in Table 1 it is clear that our phenomenological 
theory, utilizing the normal state dc resistivities as the input parameters for the microscopic 
calculation of the dielectric response, gives a good semiquantitative description of the c- 
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axis optical measurements in cuprate superconductors. Understanding the strong doping 
dependence of the data in ref. |]J requires a more complicated phenomenological model [Kj 
for t c and 7 because p c shows, in general, complex doping dependence. 

Finally, we point out that a simple phenomenological Casimir-Gorter two fluid descrip- 
tion works remarkably well for this problem. Writing the dynamical dielectric function e{u) 
as 

<«>) = 1 - ~ 4> (3) 

UJ 2 + ZCo>7 UJ 2 

where u> pn (s) is the normal (superconducting) fluid c-axis plasma frequency in the tight 
binding description obtained respectively from the normal (superconducting) carrier density, 
we can calculate (Fig. 0) the c-axis optical reflectivity. These results shown in Fig. ^| look 
qualitatively similar to our microscopic results shown in Fig. [1]. (The detailed temperature 
dependences in the two descriptions are somewhat different.) The physical picture that 
emerges from our work is thus simple: large level broadening due to strong scattering (large 
7) in the normal state associated with "incoherent" transport |10[ along the c-axis overdamps 
the normal state c-axis plasmon preventing the formation of a reflectivity edge and producing 
fairly structureless normal state reflectivity whereas the superconducting carriers, being 
unaffected by scattering effects, give rise to a sharp reflectivity edge associated with the 
well-defined Anderson-Bogoliubov c-axis plasmon which lies in the superconducting gap due 
to weak interlayer carrier hopping in the layered cuprate systems. 
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TABLES 

TABLE I. Calculated parameters t c and 7, and the corresponding c-axis plasma frequencies 
(ujpi) for various layered cuprates. 
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FIGURES 

FIG. 1. Reflectivity of LSCO as a function of the frequency (a) for fixed temperature T = OK 
(superconducting state) with various impurity scattering rate 7, (inset shows the reflectivity as a 
function of the frequency normalized by the hopping parameter t c for T = OK and 7 = A); (b) for 
fixed temperature T = 1.01T C (normal state) with various 7; (c) for fixed impurity scattering rate 
7 = A with different temperatures, and (d) for various T and 7. 

FIG. 2. (a) Reflectivity of La2- x Sr x CuO^ (x = 0.16) from the two-fluid model for different 
temperatures t = T/T c as a function of the frequency normalized by the hopping amplitude t c . (b) 
The same as in (a) for Tl 2 Ba 2 Cu0 6 (T c = 85K). 
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Table 1 



Sample 


p c (T c ) mf2-cm 


Pab (T c ) [itt-cm 


t c cm 1 


7 cm 1 


Upi (0 K) cm 1 


Upi (exp.) cm 1 


Lai.84Sro.i6Gu04 


50 [ 






140 [n] 


23 


330 


57 


55 


Tl 2 Ba 2 Cu0 6 


100 


12 




100 [0 


12 


300 


30 


< 30 [13| 


YBa 2 Cu 3 6 .6 


200 


14 


1 


100 


9 


570 


55 


60 | IS J 


Bi 2 Sr2CaCu20 8 


10 4 [ 


10 


1 


25 


0.3 


100 


5 


< 7.5 [17] 
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